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1 Introduction 


In a recent paper |l[ , all the spherically symmetric, inhomogeneous, shear-free, separable- 
metric spacetimes, that could admit matter with an equation of state p = ap where p is 
the pressure, p is the density of matter in the spacetime and ct is a constant, were reported. 
These spacetimes could be used to model a shear-free spherical star 0 or the shear-free 
universe at large 0. It is of some definite interest, and, therefore, it is the purpose of this 
paper, to investigate the Conformal Killing Vectors (CKVs) of these spacetimes. We also 
solve the corresponding geodesic equations. 


The organization of this paper is as follows. In §P we provide the salient features of 
the spacetimes of |l|. In §|3], we obtain all the CKVs for these spacetimes. In particular, 
we obtain the Killing vectors in § |3. 1| ; the Homothetic Killing vectors in jj |3.2| and the Non- 
Special Conformal Killing vectors in o[8 . 4|. There are no Special Conformal Killing vectors 


for these spacetimes as pointed out in §|373|. A summary of these results is provided in 
§[| along with the table for conformal factors for each CKV. Also provided in this section 
are the structure constants of the Lie Algebra of the CKVs. In §|5], we solve the geodesic 
equations of motion for the spacetime under consideration. In the concluding section, §j(| 
we discuss the physical applications. 


2 The Spacetime 


The metric of all the spherically symmetric, inhomogeneous, shear-free, separable-metric 
spacetimes, that admit matter with an equation of state p = ap where p is the pressure, 
p is the density of matter in the spacetime and a is a constant, is : 


ds 2 = — y 2 (r)dt 2 + 2 R 2 (t) 



dr 2 + y 2 (r)R 2 (t ) 


dd 2 + sin 2 6 deft 2 


( 1 ) 


It should be noted that the metric (P, in general, has jl|: g tt = — y 2 A 2 (t). However, 
a redefinition of the time coordinate can be made to absorb the temporal function A(t) 
without affecting any physical or geometrical properties. 

In what follows, we shall denote a time derivative by an overhead dot and derivative 
with respect to the radial coordinate by an overhead prime. Here (£, r, 9 , (j )) are the co- 
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moving coordinates. In a Lorentz frame, this metric has the Einstein tensor 


G 


tt 


G 


tr 


G rr — Ggg — G 


1 3 R 2 

2 y 2 R 2 + y 2 R 2 

y/2.R 
y 2 R 2 

1 2 R R 2 

2 y 2 R 2 ~ y^R ~ y 2 R 2 


( 2 ) 

(3) 

(4) 


Note that the G tr component of the Einstein tensor is non-vanishing for R ^ 0. The 
presence of the G tr component of the Einstein tensor implies the presence of an energy 
flux across r = constant surfaces in the spacetime of metric (|I|). This metric admits the 
following invariants: 


R = 


Z? TDdb _ 


Rn.hr.dR 


abed 


1-6 R 2 - 6RR 

Wy 2 

1 + 12i? 4 - 6 RR + 12 R 2 R 2 - AR 2 (1 - 3 RR) 

W 

3-4 R 2 + 12i? 4 - 8 RR + 12 R 2 R 2 

W 


Clearly (HD is singular iff y(r) = 0 for some r and/or R(t) = 0 for some t. 
The fluid four-velocity U a is co-moving and is given by 


U a = -5% 


1 

—< 

y 


(5) 

( 6 ) 

(7) 

( 8 ) 


The fluid four-acceleration U a = U a] f,U b is given by 

Ua = { o, —,0,0) (9) 

y 

where a semicolon denotes a covariant derivative. The spacetime under consideration has 
non-zero acceleration unless y' — 0. 

The presence of the energy flux implies^] the energy-momentum tensor of the matter 

^^While a non-vanishing radial heat flow is one possible interpretation, it is not unique. All one really 
knows is that there must be radial energy transport. This could be caused by heat, or could just as well 
be caused by the bulk motion of matter or radiation. 
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fluid to be of the form 


Tab — (p + p)U a U b + pg a b + QaUb + QbU a 


( 10 ) 


where p is the pressure, p is the density of the fluid and Q a = (0, Q , 0, 0) is the radial 
energy or the radial heat flux four-vector. Note that both shear and rotation vanish for 
the metric (]T|). 

The expansion, 0, energy or heat flux, Q, density, p and pressure, p, are given by 


0 = U a . n = 


3 R 
y R 


(ii) 


Q = 

V y 2 y'R? 


( 12 ) 


P = 


y 2 R 2 


+ 3 R 2 


(13) 


P = 


R 2 2 R 


2 y 2 R 2 y 2 R 2 y 2 R 


(14) 


From eqs. ([13]) and (|T4]), we obtain 


R = 


y 2 R 


y 2 R 2 


- (p + 3 p) 


(15) 


Therefore, the equation of state, that is, a relation between p and p, of the matter 
in the spacetime of (]lj) determines the dynamics - the behavior of the temporal function 
R(t) - of the spacetime. The temporal behavior of the heat flux is also determined by the 
temporal function R(t). 

The spatial or radial nature of the heat flux is determined primarily by the sign of 
the quantity — R/y'. The heat flux is positive, that is, heat flows from lower values of r to 

higher values of r, when y' and R have opposite signs. In the case of these two quantities 
having the same sign, the heat flux is negative, that is, heat flows from higher values of r 
to lower values of r. On the other hand, the density is, for y' > 0, a decreasing function 
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of r corresponding to a region over-dense at its center and, for y' < 0, an increasing 
function of r corresponding to a region under-dense at its center. 

For R > 0 (R < 0), the spacetime under consideration is expanding (contracting) 
since the expansion is positive (negative), 0 > 0 (0 < 0). For R > 0, the heat flux is 
positive for y' < 0 and negative for y' > 0. Therefore, when the spacetime is expanding, 
heat flows from under-dense regions to over-dense regions. For R < 0, the heat flux is 
positive for y' > 0 and negative for y' < 0. Therefore, when the spacetime is contracting, 
heat flows from over-dense regions to under-dense regions. Any two adjacent regions with 
opposite signs for heat flux are joined by heat flow caustics at locations for which y' = 0. 
Note that the quantity Q blows up at such locations. 

However, we note that the density is finite at locations y' = 0. The density is 
infinite only at locations y(r) = 0. We also observe that the metric (jl) is singular at 
heat flow caustics i.e. for y' = 0. Such locations are, however, coordinate singularities 
since the curvature invariants (j|), (|) and (]7j) are finite there. 

We also observe here that the radial function y(r) is not determined by the field 
equations. Therefore, radial attributes of matter are arbitrary , meaning, unspecified, for 
the metric (]lj). This is in the manner of concentric spheres with each sphere allowed to 
possess any value of density, for example. This is the maximal freedom compatible with 
the assumption of spherical symmetry, we may note. 

The point r = 0 will possess a locally flat neighborhood when y'\ r ~ 0 ~ 1 /x/2• 
This condition must be imposed on any y(r). Apart from this condition, the function 
y{r) is arbitrary. Other physical considerations, for example such as those arising from 
the equation of heat transfer in the spacetime, could also constrain the function y(r). 

A stellar model with heat flux can be obtained for an appropriate choice of the radial 
function y(r) [l|]. On the other hand, the metric ([[]) can also provide an inhomogeneous 
cosmological model with non-vanishing heat flow. 

This essentially completes our description of the spherically symmetric, inhomo¬ 
geneous, shear-free, separable-metric spacetimes that could admit an equation of state 
p = a p for the matter in the spacetime [Ij]. 
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3 Conformal Killing Vectors 


Conformal symmetries are of some importance in the understanding of spacetime geome¬ 
try. The conformal motions preserve the angles between vectors and, hence, the light cone 
structure of the spacetime. These symmetries thereby help solve the geodesic equations 
of motion for the spacetime under consideration. 

A Conformal Killing Vector (CKV) X satisfies 

= 2 (l>(x k ) g l3 (16) 

where <h(x fc ) is the conformal factor and g l3 is the spacetime metric tensor. There typically 
arise the following four special cases of CKVs, namely, 


1] 

Killing vectors 

Kj 

$ = 0 

2] 

Homothetic Killing vectors 

Hj 

= 0 ± $ 

3] 

Special Conformal Killing vectors 

Sj 

*ij = 0 ^ ± $ 

4] 

Non-Special Conformal Killing vectors 

Nj 



These vectors are of physical significance as they help produce first integrals. In particular, 
the Killing vectors generate the constants of motion and the homothetic Killing vectors 
scale distances by the same constant factor and, hence, preserve the null geodesic affine 
parameters. Conformal Killing vectors generate constants of motion along null geodesics. 
For details, see [[|, ||. 

The set of Conformal Killing Vectors forms a Lie Algebra G r (r < 15) with the basis 
{X;}: C Xl 9ij = 2 dy g l3 such that 


[V, X s ] = c*x k 

e 

1 

je 

= c** h 

(17) 

(16) for the metric (|lj) reduces to the system 



(M) 

^x r + x\ t 

= $ 

(18) 

(t, r ) 

-y 2 X\ r + 2 (y') 2 R 2 X r ,t 

= 0 

(19) 

C t,e ) 

-X\ e + R 2 X\ t 

= 0 

(20) 
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(t,<f>) 

-X\ 4> + R 2 sin 2 0X* t 

= 0 

(21) 

(r, r) 

|x* + yX r + XV 

= $ 

(22) 

M) 

y 2 X\ r + 2 (y'fX\ e 

= 0 

(23) 

(A0) 

y 2 sin 2 0X*V + 2{y') 2 X\ (t> 

= 0 

(24) 

(o,0) 

+ y -x r + X% 

= <f> 

(25) 

(M) 

X% + sin 2 OX^q 

= 0 

(26) 

(0» 

|x* + £x r +cot ex 9 + xV 

= <f> 

(27) 


In what follows, we obtain the Conformal Killing Vectors from the above equations case 
by case. 


3.1 Killing vectors 


As can be easily verified, the metric (|lj) admits the following three spacelike Killing vectors: 


X, 

= K r 

= (0, 0 , 0 , 1) 

(28) 

Xjj 

= Kjj 

= (0, 0, sin 0, cos (j) cot 9) 

(29) 

Xm 

= Kjjj 

= (0, 0, cos (j), — sin (p cot 9) 

(30) 


Now, for any timelike Killing vector Kjy = (K t , 0, 0, 0) orthogonal to t — constant 


hyper-surfaces, eqs. (|T9l) 


are identically satisfied. The eq. (pK|) then implies 


Xiv = K JV — ( 1 , 0 , 0 , 0 ) 


condition R = 0 


(31) 


Now, for a Killing vector Ky = (A'V), K r (r), 0, 0) that is not hyper-surface orthog¬ 
onal, eqs. (pT|), (|20|) , (pl|) , (|23|) , (p4|), ( |26|) can be satisfied. Then, eqs. 

m im Ply 


m, 


Xy =Ky 


K R, 


— kC, 0, 0 
y j 


condition R — 0 (32) 


where k is a constant and R = C. 
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The Killing vectors Kj, Kjj and Km are always spacelikc, Kjy is always timelike 
and Ky can be spacelike, null or timelike. Any spherically symmetric spacetime admits 
Kj - Km as spacelike Killing vectors. The conformal factor in each of these cases is zero. 

We note that the spacetime admitting K JV as a Killing vector is a static spacetime. 
The Killing vector K v is neither normal nor tangent to t — constant hyper-surfaces. It 
is known || that conformal motions generated by conformal Killing vectors do not, in 
general, map a fluid flow conformally. 


3.2 Homothetic Killing vectors 


In this case, <f> = constant. 


For H r — (0, H r , 0, 0), eqs. (|20|), (|2T|), ([231), (p1), ( |26| ) can be satisfied. Eqs. 

(HI), (0)’ PH), (H3)> then ’ im ply 


X VI = Hj = 


0 . bX„ 0 , 0 


(33) 


For Hjj = ( H *, 0, 0, 0), eqs. (|T1), (PD|), (PT|), (P3|), (PI), (^) can be satisfied. Eqs. 

(E3), (El)> (0) y ield 


Xyu = H n — ( B n 0, 0, 0 

XL 


condition R = 0 


(34) 


For Hjj] = (iX, 0, 0, eqs. (|T9|), (|20|), (|23|), (El), (El) can be satished. Eqs. 
(pl|), (El), (|25|), (|27| ) provide us with 


K-vm = H in — ( B in -g, 0, 0, 1 

XL 


condition R = 0 


(35) 


For Hjy = (iX, 0, H 9 , H^^j, (|T§1) and (pp fix H f . As for the angular components, there 

exist two linearly independent solutions, namely, H 6 {(j)) = sin</> and H 6 ((j)) = cos (j). We, 
therefore, obtain 


Xj X = Hjy = 


R 

B IV — , 0, sin0, cos 0 cot 6 
R 


condition R = 0 (36) 
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condition/? = 0 (37) 


( R 

By — , 0, cos 0, — sin 0 cot 6 
R 

where Bj, B II} B IU , B JV , By are constants. 

The vector /// is always spacelike, Hjj is always timelike and Hjjj, Hjy and Hy can 
be spacelike, null or timelike depending on the values of corresponding constants. The 
conformal factor in each of the above five cases is equal to the corresponding constant 
Bj, (j=i,ii,iii,iv,v). 

3.3 Special CKVs 

There are no Special Conformal Killing vectors for the metric (|Tj). The conditions $ ^ = 0, 
together with the equations ([18]) - (^7|) force $ = constant. 

3.4 Non-special CKVs 

For Non-Special Conformal Killing vectors 7 ^ 0. 

In this case we find two vectors given by 


X xr 

= Nj 

-4 t 


(38) 

Xxii 

= N n 

= i?iog 4 + ^ 

r dt d 

' 2 R dr 

(39) 


with conformal factors given by 

0 = R (40) 

0 = R logy + J 7 ^ (41) 

respectively. 

4 Summary 

It should be noted that out of the twelve CKVs given above only the Xj, Xu, Xjjj, Xyj, 
Xxi, Xxii are linearly independent. Other conformal Killing vectors reported earlier are 
obtainable from these six vectors under the stated conditions. 


In what follows, we summarize the conformal factors for these six vectors. We also 
provide in this section the Lie Brackets for these six vectors. 

The conformal factors for these CKVs are as follows: 


CKV 

Xr/Kr 

Xu/Kn 

Xrn/Krii 

$ -> 

0 

0 

0 

CKV -> 

Xyr/Hr 

Xxi/Nr 

A ' xii /Nn 

$ -»• 

Br 

R 

^logZ/ + /|| 


The above CKVsp] form the basis of the Lie Algebra of CKVs for the metric (]!]). 
The Lie Algebra of CKVs is therefore six dimensional, that is, G&. The Lie brackets and 
corresponding structure constants are easily obtained as follows: 

[K r , Kj T \ = Kjjj 

[Ku,K in } = Kj 
[K r jj,K r \ = K n 
\Hj,N n ] = N r 

[Nj,N n ] = l -Hj 

We note that the three Killing vectors - Kj, Kn, Km - form a sub-algebra corre¬ 
sponding to the rotational symmetry of the space-time. On the other hand, the remaining 
three Conformal Killing Vectors - Hi , Nj, Nji - form another sub-algebra isomorphic 
to the Lie Algebra corresponding to the group E + ( 2) ||. This is realized by setting 

Ui = Hr, IJ 2 = — \/2 i Nr , V 3 = \[2iNn. This is an interesting aspect of the space-time 
of the metric ([!]). 

2 Note: The above results have been confirmed using the differential equation solver within PROGRAM 
LIE 0. 
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5 Geodesic equations of motion 


The Lagrangian for the metric ([!]) i 


is 


2 C = — y 2 1 2 + R 2 2 (y') 2 r 2 + y 2 6 2 + y 2 sin 2 9 12 


(42) 


where an overhead ~ denotes a derivative with respect to the affine parameter s along 
the geodesic. Further, C takes values 0 for null, +1 for spacelike and — 1 for timelike 
geodesics. 

The geodesic equations of motion are obtained from the Euler-Lagrange equations: 


d [3(2C)\ <9(2 C) 


ds \ dx a 


dx° 


(43) 


In what follows, we will use (^2]) to reduce the terms after rearrangement. 
The r-equation is obtained as: 


d_ 

ds 


( R 2 yy’r ) = C 


and can be easily integrated to obtain 


r = 


C s + ki 
Bfyy' 


(44) 


(45) 


The 0-equation is 


and it implies 


— (-R 2 y 2 sin 2 6^ =0 


ko 


R 2 y 2 sin 2 6 


where k is a constant of integration. 
The 0-equation is obtained as 


(46) 


(47) 


( i? 2 y 2 0 ) = R 2 y 2 sin 0 cos 0 0 2 


ds 


(48) 
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On using ( 0 ) and rearranging terms, we obtain 


d 

ds 


(r 2 y 2 &') + k\ cot 2 6 


= 0 


(49) 


which implies 

(i ? 2 y 2 d^j + fc 2 cot 2 9 — k, 
where k 3 is an integration constant. Hence, 


0 = \j h ~ kl cot 2 6 

R z y z v 


(50) 


The t-equation is obtainable as 


ds 


( Ry 2 ~t ) 


= -2CR 


(51) 


where we have used (|42|) to reduce terms. 


It is difficult to integrate this equation directly. However, the solution is obtainable 
from (|42|) using (|5|), (|7]), (H) as 


t = 


9 C 1 

+ —[2 ( jCs + s ) 2 + k 3 + kf\ 


y 2 y A R 2 


(52) 


Therefore, the geodesic of the metric ([!]) has the tangent vector 


rpa 


( i, f- ft i ) 


(53) 


with t, r, 9, cj) as obtained above. 


Of course, for null geodesics, we obtain by direct integration of (|5l|): 


Ry 2 t = 5 


(54) 


where 6 is a constant. Therefore, in general, the null geodesic of the metric (Ij) has the 
tangent vector 


n = 


y 2 R 2 


(sR, ±\Jk 3 - k 2 cot 2 6, - 7 %^) 

y y ' v sm 9 j 


(55) 
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with the different constants satisfying a relation obtainable from n a n a = 0. In particular, 
for the radial null geodesic we have 

5 = ±V2h (56) 

Along a null geodesic x a (s ) with tangent vector n a = dx a /ds , the following holds: 

n a - b n b = 0 = n a n a (57) 

Any conformal Killing vector, X, generates along a null geodesic a constant of the motion: 

y-(X«n) = X a . b n a n b = 2 <5>g ab n a n b = 0 (58) 

ds 

Each geodesic is specified by six parameters: four to determine a point on it and two 
to determine the direction of its tangent. Thus, the general solution of the null geodesic 
equation for the metric ([I]) depends on six functionally independent constants of null 
geodesic motion. These constants of motion are obtainable from the six conformal Killing 
vectors and the above tangent to the null geodesic trivially. 

6 Discussion 

Most of the known solutions of the Einstein field equations admit some nontrivial isom¬ 
etry group. Further, there are many solutions with conformal symmetry, but in most of 
these the symmetry is actually homothetic. Much physical insight on astrophysical and 
cosmological questions has been obtained from the study of such solutions. Of course, 
with the assumption of the spherical symmetry, there is the rotational isometry group 
associated with the spacetime. The question then is of obtaining further symmetries of 
the spherically symmetric spacetime. In this respect, our spacetime here is interesting 
since it has only one homothetic symmetry but two conformal symmetries. A homothetic 
Killing vector leads to self-similarity while scaling all the distances by the same constant 
factor. The conformal symmetries provide the constants of geodesic motion along null 
trajectories of massless particles. 

In this paper, we obtained all the conformal Killing vectors for the metric (|l[) and 
their Lie Algebra. We showed that, for the space-time of the metric ([!]), there exists a sub¬ 
algebra isomorphic to the Lie Algebra of E + { 2). This Lie Algebra has an associated Local 
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Lie Transformation Group E + ( 2) that is a group of distance preserving transformations 
in the plane i?, 2 . 

We also obtained explicitly the tangent vector to the geodesics of the metric ([[]) . In 
particular, for the trajectories of massless particles or the null geodesics, the constants 
of motion can be obtained from the CKVs presented. The results obtained here are of 
importance for the study of astrophysical problems involving stellar models or for the 
study of inhomogeneous cosmological models based on the metric (/!]). 
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